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ABSTRACT 

We present RXTE observations of two recently identified massive X-ray bina- 
ries. RX J0440.9+4431/BSD 24-491 and RX J1037.5-564/LS 1698 are confirmed as 
accreting Be/X-ray systems following the discovery of X-ray pulsations, with barycen- 
tric pulse periods of 202.5±0.5 s and 860±2 s respectively. The X-ray spectral analysis 
shows that the energy spectra of the pulsars can be represented by a power-law, mod- 
ified at low energy by an absorption component and at high energy by a cut-off. Very 
weak Fe lines may be present. Both sources appear to display a low cut-off energy 
when compared to typical X-ray pulsars, low X-ray variability (factor of <10), and no 
dependence of the X-ray spectrum with energy. Given the similarity of these X-ray 
properties with those of the other persistent BeXRB pulsars, 4U0352+309/X Per and 
RX J0146.9+6121/LS 1 +61 235, we suggest that RX J0440.9+4431/BSD 24-491 and 
RX J1037.5-564/LS 1698 are also members of this subclass. 
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1 INTRODUCTION 

Motch et al. (1997) (hereafter M97) reported the discov- 
ery of several new high mass X-ray binaries. These sys- 
tems were found during the ROSAT galactic plane survey by 
cross-correlating the position of low-latitude X-ray sources 
(|6| < 20°) with SIMBAD OB star catalogues. The search 
was restricted to stars earlier than B6 and X-ray luminosi- 
ties > 10 31 erg s _1 . When the X-ray data were combined 
with optical observations, five sources remained candidates 
for new accreting binary systems. These sources are BSD 
24-491 (RX J0440. 9+4431), LS 992 (RX J0812. 4-3114), LS 
1698 (RX J1037.5-5647, also probably 4U1036-56 accord- 
ing to M97), LS 5039 (RX J1826. 2-1450) and LS I +61 235 
(RX J0146.9+6121). In this paper we present RXTE obser- 
vations of BSD 24-491 and LS 1698. Hereafter, we will refer 
to sources by the names of the optical companion. 

Be/X-ray binaries (BeXRBs) and supergiant systems 
represent the general class of High Mass X-ray Binaries. 
These are systems consisting of an early type star, the ex- 
act size forming part of the classification of the system. In 
addition to the normal stellar component, a compact object 
is also present in the system, and accretion onto it is the 
principal source of X-ray emission in the system. Follow- 
up optical spectroscopic observations carried out by M97 
showed Ha in emission in both, BSD 24-491 and LS 1698, 



hence classing them as Be stars. The spectral type for both 
sources is BOV-IIIe (M97). 

A Be star is an early type luminosity class III-V star, 
which at some time has shown emission in the Balmer series 
lines. This emission, as well as the characteristic infrared ex- 
cess, is attributed to the presence of circumstellar material, 
most likely forming a disc around the equator of the Be star. 

BeXRBs were first identified by the presence of X-ray 
flares, and so were named X-ray transients. Type I flares 
were observed to occur regularly, with a fixed period, and it 
was hypothesised that the systems consist of a neutron star 
in an eccentric orbit around the Be star. This model predicts 
that an X-ray flare will be observed during the time of the 
neutron star's periastron passage, and this explains the pe- 
riodicity observed in the timing of the flares. The size of the 
flare will be dictated by the amount of material available for 
accretion and by the magnetic field of the neutron star. The 
flare occurs at periastron when the neutron star impinges 
on the region of denser circumstellar material surrounding 
the Be star. BeXRBs also show giant X-ray outbursts (Type 
II events) in which the X-ray luminosity increases by a fac- 
tor 100-1000 above the quiescent level. The X-ray luminos- 
ity during the outburst may reach 10 38 erg s _1 , close to the 
Eddington luminosity. These outbursts do not correlate with 
the orbital phase, but occur in an unpredictable way and are 
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thought to be related to mass ejections from the Be star's 
circumstellar envelope. 

It has been proposed that a subclass of BeXRBs exists 
which are characterised by persistent, weak X-ray emission, 
and which do not display the type I or II outburst behaviour. 
This subclass is currently represented by X Persei and LS I 
+61 235 (Haberl, Angelini, Motch and White, 1998). These 
systems have relatively low X-ray luminosities (~10 34 erg 
s _1 ) and long spin periods (837 s X Persei and 1412 s LS I 
+61 235). To date, these two represent the only examples of 
persisent BeXRBs amongst the ~40 known BeXRB systems. 
For a recent summary of the observed properties of BeXRBs, 
see Negueruela (1998). 



2 OBSERVATIONS 

2.1 X-ray observations 

The sources were observed with the Proportional Counter 
Array (PCA) onboard the Rossi X-ray Timing Explorer 
(RXTE) in 1998 February. Table 1 shows the journal of the 
X-ray observations. The total on-source time was 20 ks for 
each target. The PCA covers the lower part of the energy 
range, 2-60 keV, and consists of five identical coaligned gas- 
filled proportional units (PCU) , providing a total collecting 
area of ~ 6500 cm 2 , an energy resolution of < 18 % at 6 
keV and a maximum time resolution of lfis. For a more 
comprehensive description of the RXTE PCA see Jahoda et 
al. (1996). 

Good time intervals were defined by removing data 
taken at low Earth elevation angle (< 8°) and during times 
of high particle background. An offset of only 0.02° between 
the source position and the pointing of the satellite was al- 
lowed, to ensure that any possible short stretch of slew data 
at the beginning and/or end of the observation was removed. 
This screening criteria allowed us to divide the observations 
up into continuous sections of clean data, on which the X- 
ray analyses were carried out. The main objective of these 
observations was to search for pulsations in the X-ray flux in 
an attempt to confirm the accreting nature of the sources. 



3 TIMING ANALYSIS 

In this section we present the results of the X-ray timing 
analysis. For each source lightcurves, pulse profiles and hard- 
ness ratios are given and the best value of the spin period 
determined. 



3.1 BSD24-491 

Fig hi shows an enlarged section of the lightcurve at three 
different energy ranges: 3-6, 6-10, 10-20 keV. Pulsations are 
clearly seen at the three energy ranges considered and they 
show up in the Power Density Spectrum (PDS) as a peak 
at ~ 0.005 Hz, which corresponds to a spin period of ~ 
200 s. The X-ray intensity of the source, after background 
subtraction, was 10.7+0.1 PCA c s _1 in the energy range 3- 
30 keV. Other than pulsations, no other type of variability 
was seen during the ~ 6 hour interval of our observation. 
The extrapolated luminosity in the energy range 0.1-2.4 keV 
is 1.4 x 10 33 erg s _1 , compared to the maximum value of 



Figure 1. X-ray pulsations in BSD24-491 at three different en- 
ergies, with each bin representing 16s. The starting time is JD 
2,450,843.398. Also shown is the variation of the HR 10-20/6-10 
keV with time 



6.0 x 10 33 erg s" 1 given by M97 during a ROSAT pointed 
observation, indicating some long term variability. In fact 
M97 reported a variation in the count rate of a factor of 2.5 
over their observations. 

In order to determine the pulsation period, a period 
search by epoch folding was performed near to the period 
expected from the FFT power spectrum. A solar barycentric 
pulse period of 202 s was obtained from the resulting pulse- 
folding periodogram. The pulse profile obtained from this 
period was used as a template. Then, both the observations 
and the template were divided into 4 roughly equal time 
intervals. The difference between the actual pulse period at 
the epoch of observation and the period used to fold the 
data was determined by cross-correlating the original light 
curve and the template. The P actua i - Ptemplate shifts were 
fitted to a linear function to obtain an improved period. A 
new template was derived from this period and the process 
repeated. The derived pulse period was 202.5+0.5 s. The 
error represents the scatter of the points about the best-fit 
straight line. The accuracy in the determination of the spin 
period is limited by the length of the time series with respect 
to the length of the period. 

Once the pulse period is known, the pulse profile can be 
derived by folding the lightcurve modulo the best fit value 
of the pulse period. Such a pulse profile is shown in Fig @ 
at the same three energy ranges in which the lightcurve was 
separated. No depence on energy is seen. The pulse shape is 
nearly sinusoidal, as might be expected from the absence of 
second or higher harmonics in the PDS. The amplitude of 
the modulation is ~ 75% in the 3-30 keV energy range. 
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Table 1. Journal of the X-ray observations 



Name of 


Date 


Start time 


Stop time 


Luminosity" ' b 


Spin Period 


HR C 


Source 








(erg s" 1 ) 


(s) 




BSD 24-491 


30/01/98 
01/02/98 


06:24:25 
05:26:38 


11:19:14 
11:07:14 


3.0X10 34 


202.5±0.5 


-0.56 


LS 1698 


04/02/98 
05/02/98 

05/02/98 


16:25:14 
08:40:50 
10:03:38 


19:49:14 
09:20:14 
16:35:14 


4.5 X10 35 


860±2 


-0.40 



a: in the energy range 3-30 keV 

b: The assumed distances are (from M97) BSD 24-491: 3.2 kpc, LS 1698: 5.0 kpc 

c: HR=(10-30)-(3-10)/(3-30) 



Figure 2. Pulse profiles of BSD 24-491 at three different energy 
bands 



3.2 LS 1698 

M97 reported a maximum X-ray luminosity L x (0.1-2.4 keV) 
« 1.1 x 10 34 erg s _1 for LS 1698. The average 3-30 keV 
RXTE PCA count rate, after background subtraction, was 
45.8±0.1 count s _1 , without significant changes in intensity. 
The PDS averaged over two intervals of 3600 s revealed the 
presence of one peak at v ~0.0011 Hz, corresponding to a 
coherent modulation with a period of ~900 s. 

In order to determine a more accurate value of the 
pulse period we rebinned the lightcurve into 32 s bins 
spanning a total on-source time of ~ 23500s (including 
gaps due to Earth occultations or passage through the 
South Atlantic Anomaly) . Times were corrected to the Solar 
System bary centre, and periodicities searched for in these 
lightcurves using the epoch folding technique. Data were 
folded over a range of trial periods and the \ 2 0I the folded 
lightcurve calculated. The period adopted is the one which 
shows a maximum in the x 2 versus period diagram. The 
period providing the highest \ 2 was found to be 860±2 
s. Although the observations do not cover a long enough 
base line to further constrain this period, e.g. by applying 
the cross-correlation method explained above, pulsations arc 
distinctly seen in the X-ray lightcurve (Fig H). Fig W dis- 



plays the 3-6, 6-10 and 10-20 keV folded lightcurves of LS 
1698. These pulse profiles are highly structured with a broad 
peak covering phases 0.1-0.7. The pulse fraction, PF=(I max - 
Imin) / 0-max+lmin) is ~53±2% and remains the same for the 
three energy ranges considered. I max and I min are the inten- 
sity at pulse peak and pulse dip, respectively. 

M97 noted that the source 4U1036-56/3A1036-565 may 
be the same as RX J1037.5-564/LS 1698, based on the po- 
sitional coincedence of the ROSAT error circle with those of 
Uhuru and Ariel V. However, despite observations with sev- 
eral satellites (Uhuru, Ariel V, OSO 7, EXOSAT, ROSAT), 
no pulsations have previously been detected from 4U1036- 
56 (although this may simply be attributable to the lower 
sensitivity of previous observations rather than a real indica- 
tion that no pulsations are present). The reported detection 
of a flare in the Ariel V observations (Markert et al. 1979) 
only shows a slight increase (x2.4) over the Uhuru flux in 
the same energy region (2-10 keV), and the ROSAT obser- 
vations were about an order of magnitude lower than this. 
The Uhuru and Ariel luminosities are consistent with those 
we obtain from our RXTE observations, ~3 x 10 35 erg s _1 
compared to our value of ~2 x 10 35 erg s _1 (2-10 keV and 
assuming a distance of ~5 kpc). It would thus appear that 
4U 1036-56 displays a low level of X-ray variability, which 
might support the identification with LS 1698, but we can- 
not confirm this. 



3.3 Hardness ratios 

LS 1698 is the hardest of the two sources, as indicated by 
a hardness ratio (HR) of -0.40, defined as (10 — 30) — (3 — 
10)/(3 - 30), compared to -0.56 displayed by BSD 24-491 
(see Table 1). This fact also becomes apparent in Fig 0, 
which shows tha hardness ratios in the photon-energy bands 
10-20 and 3-6 keV as a function of the summed count rate 
of the two bands. It is worth pointing out a few interest- 
ing results. Firstly, the HR of BSD 24-491 and LS 1698 re- 
mains basically constant at ~0.4 and ~0.9 respectively. This 
contrasts with the behaviour of the 'typical' BeXRB pulsar 
transient system RX J0812-3114/LS992 shown for compar- 
ison, where a clear dependence on intensity is seen (Reig & 
Roche 1999). Secondly, as stated above, the X-ray emission 
from LS 1698 is much harder than that of BSD 24-491. 
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Figure 3. X-ray lightcurve of LS 1698 showing pulsations at 
different energy bands. Each bin represents 32 s. Time is JD 
2,450,849.420. Also shown is the hardness ratio 10-20/6-10 keV 
as a function of time. Note that the HR follows the pulses. 



Figure 5. Hardness ratio 10-20/3-6 keV as a function of intensity 
of the persistent BeXRBs BSD 24-491 and LS 1698. Shown also 
for comparison is the HR variation of a transient BeXRB (from 
Reig & Roche, 1999) 



- broad-band emission, not dominated by sharp spectral 
features 

- a continuum described in terms of a power-law with low 
energy absorption and a cut-off above 10-20 keV 

- the X-ray flux is mostly emitted in the energy range 
2-20 keV with a rapid falloff above 20 keV 

- an iron line at ~6.4 keV due to fluorescent reprocess- 
ing by cold circumstellar material is present in a number of 
systems 

In order to investigate the X-ray energy emission of our 
sources, fits using a variety of models were performed. Table 
2 gives the results of the spectral analysis. We have included 
only the values of the models which gave acceptable fits, i.e. 



reduced \ < 2. 



Figure 4. Pulse profiles of LS 1698 at 3-6, 6-10 and 10-20 keV. 



4 SPECTRAL ANALYSIS 

Although each source adds its peculiarities, the X-ray pho- 
ton energy spectra of accreting X-ray pulsars share the fol- 
lowing general characteristics: 



4.1 BSD 24-491 

The analysis of the X-ray spectrum of BSD 24-491 is compli- 
cated by the low count rate. To improve the signal to noise 
ratio we used only data from the top layer anodes of the 
four PCUs that were on during the entire observation. Also, 
models were fitted to the energy range 3-20 keV instead of 
the 3-30 keV used for LS 1698. In spite of this, the errors are 
too large to meaningfully distinguish between the different 
spectral models listed in Table 2. An F-test shows that the 
inclusion of a blackbody component to the power-law model 
is not justified by the reduction of the % 2 - F° r the purpose of 
the subsequent discussion, and in order to compare spectral 
parameters with the other sources, we will consider only the 
results from the power-law plus cut-off model. An emission 
feature at around 6.2 keV shows up in the residuals (Fig 
ph but the fit does not allow us to constrain its parameters 
precisely. We set an upper limit for the equivalent width of 
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the iron line of ~ 100 eV. The line centre remained at the 
same value of 6.2±0.2 keV despite fixing the line width to 
a ~ 0.1, 0.2 and 0.3 keV, indicating that the line width is 
narrower than the spectral resolution of the PC A (18% at 
6 keV) . BSD 24-491 also shows a very low cut-off energy of 
~1.9 keV. This value is outside the detector response range 
and is thus not well constrained. However, if the iron line 
is included in the fit the cut-off energy increases to ~4.5 
keV. In the energy range 3-30 keV, L x ~3.0 x 10 34 erg s _1 , 
assuming a distance of 3.2 kpc. 
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4.2 LS 1698 

The best fit to the X-ray spectrum of LS 1698 is obtained us- 
ing the power-law plus high energy cut-off model (x?— 1-23 
for 55 dof). The inclusion of a Gaussian component, repre- 
senting an iron line, improved the fitting by reducing \r to 
1.00 for 52 dof. However, an F-test shows that the proba- 
bility of this happening by chance is higher than 10%. The 
best-fit line centre energy and width are 6.5±0.2 keV and a 
~0.2 respectively. If the Gaussian component is taken into 
account, the cut-off energy is 4.7 keV, much lower than the 
values of 10-20 keV typically found in X-ray pulsars. This 
may represent a signature of the persistent BeXRB systems 
(see below) . An upper limit for the equivalent width of the 
iron line is estimated to be ~65 eV. The photon index lies 
in the interval a=1.0-1.2 depending on whether the iron 
line is included in the fit or not. LS 1698 has E(B-V)^0.75 
(M97), implying N H ~0.5 x 10 22 cm~ 2 (Ryter, Ceasarsky 
& Audouze, 1975). However, the absorption implied from 
the spectral fitting is higher, ~4.5 x 10 22 cm -2 . Thus, some 
amount of matter must be located close to the system, which 
would justify the presence of the iron line. The X-ray lumi- 
nosity is ~4.5 x 10 35 erg s _1 in the energy range 3-30 keV 
and assuming a distance of 5 kpc. 



5 DISCUSSION 

Pulsations have previously been found in both high and 
low mass X-ray binaries. However, of the 38 X-ray pulsars 
currently known, only 5 have been identified as LMXRBs, 
whereas 33 have an OB-type star as the optical counter- 
part. For the confirmed BeXRB systems, around 70% are 
found to be pulsators, rising to 75% if suspected BeXRBs 
are included. 

The periods of X-ray pulsars are distributed over a fac- 
tor ~10 5 , from 2.5 ms (SAX J1808-369, the bursting mil- 
lisecond X-ray pulsar) to 1412 s (LSI +61 235, a persistent 
BeXRB), with no evidence for clustering at any particular 
period. The range for BeXRBs is almost as extensive, the 
fastest being 69 ms (A 0535-668). 

We have carried out X-ray timing and spectral analyses 
of two of the four newly discovered high mass X-ray binaries 
by M97. We have detected pulsations in both BSD 24-491 
(period 202.5 s) and LS 1698 (~860 s). In addition to slow 
pulsars, LS 1698 and BSD 24-491 seem to share further sim- 
ilarities in their X-ray variability characteristics. They have 
both been detected every time that they have been observed 
(during the ROSAT all-sky survey and subsequent follow- 
up pointed observations using ROSAT (M97) and here with 



Figure 6. PCA spectrum of BSD 24-491. The continuum was fit 
to a power-law plus a high energy cut-off. If no Gaussian compo- 
nent is added to the fit, an emision feature at sy 6.2 keV shows 
up in the residuals 



RXTE) and they do not appear to show short-term variabil- 
ity (of time scales of hours or days). 

5.1 Is there a class of persistent Be/X-ray 
binaries? 

Until recently, X Per/4U0352+309 was unique among mas- 
sive X-ray binary systems in that it was the only persistent 
BeXRB. The optical component is a BOVe star at a distance 
of ~700 pc (Lyubimkov et al. 1997). Persistent high mass 
X-ray binaries do exist but they harbour an evolved, i.e. 
supergiant, companion. In these systems the X-rays are at- 
tributed to constant accretion from the strong stellar wind 
of the primary. 

In the case of X Per, the weak but persistent emission 
is most likely a result of the fact that the neutron star is 
always a long way from the Be star and never enters the 
denser inner regions of the circumstellar disc. An orbital pe- 
riod of 580 days has been proposed (Hutchings et al. 1974), 
but never subsequently confirmed by later studies. However, 
there appear to be characteristic timescales of photomet- 
ric variability ranging from hundreds to thousands of days 
(Roche et al. 1997), but nothing that appears to be orbital 
in nature. In most respects, X Per behaves like a classical Be 
star, presumably because the separation of the components 
means that the B0V star and it's disc are unaffected by the 
distant neutron star. The neutron star can be affected by 
episodes of mass ejection from the circumstellar disc (e.g. 
Roche et al. 1993, 1997), but is usually only accreting from 
the Be star wind. 

Recently, Haberl et al. (1998) and Haberl, Angelini 
& Motch (1998) have suggested that LS I +61 235/RX 
J0146.9+6121 may be another persistent BeXRB. This 
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Table 2. Spectral fits results. Uncertainties are given 90% confidence for one parameter of interest. Spectra were fitted in the energy 
range 3-30 keV for LS 1698 and 3-20 keV for BSD 24-491. 





BSD 24-491 


LS 1698 


Power-law 






a 


2.33+0.07 




N H (10 22 atoms cm" 2 ) 


6.2±0.7 




X?(dof) 


1.14(43) 




Power-law & blackbody 






a 




1.84+0.06 


kT (keV) 




2.9+0.2 


R (km) 




0.10+0.01 


N H (10 22 atoms cm" 2 ) 




8.2±0.5 


xHdof) 




1.66(54) 


Two blackbody 






kTi (keV) 


1.30±0.06 


1.55±0.04 


Rl (km) 


0.22±0.01 


0.49+0.02 


kT 2 (keV) 


3.8±0.6 


4.48+0.12 


R 2 (km) 


0.03±0.01 


0.085±0.004 


N H (10 22 atoms cm" 2 ) 


- 


2.1±0.3 


xHdof) 


0.78(43) 


1.34(55) 


Cut-off power-law * 






a 


1.5±0.4 


1.20±0.08[1.02±0.07] 


Ecut 


~ 1.9[~4.5] 


6.2±0.4[4.7±0.4] 


N H (10 22 atoms cm" 2 ) 


4.1±1.0[~1.9] 


5.5±0.5[4.5±0.4] 


Xr(dof) 


.94(41) [0.56(40)] 


1.23(55)[0.98(53)] 



* The values in 



represent the best-fit parameters when an iron line is included. For BSD 24—491, E c 
keV, whereas for LS 1698 E c =6.5±0.2 keV and a sj 0.2 keV 



=6.2±0.2 keV with fixed cr=0.1 



Table 3. Comparison of the X-ray properties of persistent and transient BcXRBs 



X-ray 
name 



Optical 
name 



Spectral 
type 



Spin 
period (s) 



Spectral 
index 



E cu t 
(keV) 



T 6 

-'-'x.max 

x 10 36 



X-ray EW(Fe) 

variability (cV) 



4U0352+309 
RX J0146. 9+6121 
RX J1037.5-5647 
RX J0440.9+4431 



X Per BOVe 

LS I +61 235 BlVc 

LS 1698 BOV-IIIe 

BSD 24-491 BOV-IIIe 



837 


0.8+0.2 


2.2+0.4 


0.03 


5 


< 100 


1404 


1.2+0.2 


4.0+0.7 


0.5 


20 


< 90 


860 


1.0+0.1 


4.7+0.4 


0.5 


20 


< 65 


202 


1.5+0.4 


~4.5 


0.03 


6 


< 100 



EXO 2030+375 
4U0115+63 



V635 Cas 



Be 
09e 



41.7 
3.6 



1.3-1.8+0.01° 
0.36+0.01 



5-19+0.5° 
7.5+0.2 



100 
30 



1000 
200 



200 
110 



a dependent on luminosity; higher values correspond to higher luminosity 
b energy range 1-20 keV except that of LS 1+61 235 which is 0.5-10 keV 

C l^x,max/t^x,min 



source is the pulsar with the longest spin period (~ 1400 
s) . With two such systems one can look for more similarities 
in an attempt to establish the existence of a distinct subclass 
of BeXRBs. Comparing the X-ray observational characteris- 
tics of X Per and LS I +61 235 we find that they both share 
the following properties: 

- Long pulse periods, 837s for X Per and 1404s for LS I 
+61 235 

- Persistent, low luminosity (<10 34 ~ 35 erg s" 1 ) X-ray 
emission 

- Low cut-off energy derived from the spectral fitting, 
compared to the typical value of 10-20 keV found in transient 
X-ray pulsars - 2.2 keV for X Per (Schlegel et al. 1993) and 
4 keV for LS I +61 235 (Haberl, Angelini & Motch 1998) 



- Absent or very weak iron line at 6.4 keV, indicative of 
only small amounts of material in the vicinity of the neutron 
star 

- Low X-ray variability. Flat lightcurves with rare and 
unpredictable increases in flux by a factor of < ~10 

- No dependence of the X-ray spectrum on intensity 

If we use these general characteristics as defining a sub- 
class of BeXRBs, we find that the newly discovered systems 
LS 1698 and BSD 24-491 are also potential members. Both 
have been detected every time they have been observed, in- 
dicating persistent X-ray emission, and their X-ray flux is 
relatively low compared to that of other accreting pulsars 
when they are active. Furthermore, both have long spin pe- 
riods (-860 s for LS 1698, ~202.5 s for BSD 24-491) and 
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Figure 7. PCA spectrum of LS 1698 and residuals. The best-fit 
power-law with high energy cut-off is plotted as a straight line 



combination of little material being accreted plus slow ro- 
tation plus possibly a relatively weak magnetic field makes 
these systems weak but persistent pulsing sources. 

Due to the low accretion rate outside of these rare events 
(e.g. X Per has undergone only one observed episode of X- 
ray flaring activity in almost 30 years, associated with major 
changes in the circumstellar environment, Roche et al. 1993 
and 1997; LS I +61 235 only 2 flares in 13 years, Haberl, 
Angelini & Motch 1998), the neutron star may not be spin- 
ning at the equilibrium period, and will thus not follow the 
relationship predicted from the Corbet diagram of P 3p i n v 
P orbit (Corbet 1986; Waters & van Kerkwijk 1989). How- 
ever, a long orbital period is expected, to account for the 
presumed distance of the neutron star from the envelope. 

Of these persistent systems, only X Per has been the 
subject of detailed, long-term study. In this case, whilst 
many of the physical characteristics of the system have 
been elucidated (e.g. stellar parameters by Lyubimkov et 
al. (1997), disc parameters by Telting et al. (1998), photo- 
metric and polarimetric behaviour by Roche et al. (1997) 
etc.), there are still doubts about fundamental parameters 
such as the orbital period. None of the proposed members 
of this subclass has an orbital period determination, and so 
their position on the Corbet diagram cannot be ascertained. 
Future work should therefore focus on a search for orbital 
variations in these objects. 



low cut-off energies of 6.2 (4.7) and 1.9 (4.5) keV respec- 
tively - the values in brackets are obtained when an iron 
line is included in the fit, but this is not statistically re- 
quired, again as in the cases of X Per and LS I +61 235. We 
have also shown in Sect 3.3 that the HRs did not change sig- 



nificantly throughout the observations, indicating that the 
X-ray spectrum is insensitive to changes in the count rate. 
This absence of spectral changes with intensity was also seen 
in LS I +61 235 by Haberl, Angellini & Motch (1998). In 
contrast, transient BeXRBs show, in general, a positive cor- 
relation between HR and intensity, as seen in LS 992 (Reig 
& Roche 1999) or EXO 2030+375 (Reig & Coe 1998). It is 
worth noting that displaying just one of this characteristics 
is not enough to qualify as a member of the group, but the 
combination is suggestive of a real association. For example, 
the X-ray spectrum of LS 992 also presents a low cut-off 
energy (4.7 keV), but it does not share any of the other 
X-ray properties (Reig & Roche, 1999). The X-ray proper- 
ties of these persistent BeXRBs and some examples of the 
transient BeXRBs are compared in Table 3. 

It therefore seems that we now have a growing subclass 
of BeXRBs, characterised by persistent, low-luminosity X- 
ray emission and slowly rotating pulsars. A possible model 
for these systems is that of a neutron star orbiting a Be 
star in a relatively wide orbit, accreting material from only 
the low density outer regions of the circumstellar envelope. 
Sporadic ejection of the disc itself, however, will result in the 
ejected material and the neutron star interacting, producing 
a flare in the X-ray lightcurve. The increase in flux during 
these flares is typically less than a factor ~10 (Roche et al. 
1993), although further work needs to be done on character- 
ising the range of X-ray variabilities observed in the whole 
range of BeXRBs. Note that the expelled material can reach 
the neutron star at any orbital phase (not necessarily at pe- 
riastron) and that flares need not be periodic. Thus, the 



6 CONCLUSIONS 

We have discovered two new Be/X-ray pulsars which ap- 
pear to display X-ray characteristics in common with the 
proposed subclass of persistent BeXRBs currently consist- 
ing of X Per and LS I +61 235. If the association of LS 
1698 and BD 24-491 with these other objects is correct, the 
sample of persistent BeXRBs has doubled, and further ob- 
servations of all these sources is encouraged to study the 
properties of this subclass. Currently, only X Per has been 
studied extensively for any duration, but the existence of 
similar systems provides us with probes of the behaviour of 
neutron stars in low density accretion regimes. A possible 
model for these systems is that of a neutron star orbiting 
a Be star in a relatively wide orbit, accreting material from 
only the low density outer regions of the circumstellar en- 
velope. Small increases in X-ray luminosity (e.g. of order 
10 times the quiescent flux) may result from mass ejection 
episodes where the neutron star is temporarily in a region 
of increased matter density. 
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